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Abstract 

Understanding traffic statics and dynamics in urban networks is critical to develop effective 
control and management strategies. In this paper, we provide a novel approach to study the 
traffic statics and dynamics in a signalized double-ring network, which can provide insights 
into the operation of more general signalized traffic networks. Under the framework of the link 
queue model (LQM) and the assumption of a triangular traffic flow fundamental diagram, the 
signalized double-ring network is studied as a switched affine system. Due to periodic signal 
regulations, periodic density evolution orbits are formed and defined as stationary states. A 
Poincare map approach is introduced to analyze the properties of such stationary states. With 
short cycle lengths, closed-form Poincare maps are derived. Stationary states and their stability 
properties are obtained by finding and analyzing the fixed points on the Poincare maps. It 
is found that a stationary state can be asymptotically stable, Lyapunov stable, or unstable. 

The impacts of retaining ratios and initial densities on the macroscopic fundamental diagrams 
(MFDs) and the gridlock times are analyzed. Multivaluedness and gridlock phenomena as well 
as the unstable branch with non-zero average network flow-rates are observed on the MFDs. 

With long cycle lengths, fixed points on the Poincare maps are solved numerically, and the 
obtained stationary states and the MFDs are very similar to those with short cycle lengths. 

Compared with earlier studies, this paper provides an analytical framework that can be 
used to provide complete and closed-form solutions to the statics and dynamics of double¬ 
ring networks. This can lead to a better understanding of how the combination of signalized 
intersections and turning maneuvers is expected to impact network properties, like the MFD. 
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1 Introduction 


The rapid increase in travel demands has led to severe congestion problems in urban transporta¬ 
tion networks ( [Papageorgiou et aQ|2003| >. This congestion has many significant negative impacts 
to society, including lost productivity, wasted fuel, and the creation of environmental and health 
problems ( |Schrank et al412012| >. One recommended strategy to alleviate urban traffic congestion 
is improving signal timings at intersections where the current technology is deficient ( [Sorensen 
et al.[ |2008| >. Many control strategies have been proposed since one of the earliest attempts in 
( |Webster[ |1958| >, but most are developed for under-saturated conditions and are not applicable to 
over-saturated conditions ( [Papageorgiou et aL| |2003[ ). In order to develop effective control and 
management strategies for urban networks, a fundamental challenge is to understand the static and 
dynamic properties of traffic in signalized networks under both under- and over-saturated condi¬ 
tions and a variety of other factors, e.g., demand patterns, signal settings, and route choice behav¬ 
iors. 

To evaluate the impact of changes (e.g., re-timing the signals) on urban networks, one practi- 


cal approach is via the macroscopic fundamental diagram (MFD)( Daganzo and Geroliminis, 

2008 

). 

The MFD is a network-level flow-density relation, which was first proposed by ( 

Godfrey, 

1969 

). 

Recent studies have shown that MFD exists in urban networks (Geroliminis and Daganzo[[2008) 

and is related to many factors, such as spatial variability in congestion and density distributions 

(Buisson and Ladier 

, 2009[ Ji et al.[ 2010; Geroliminis and Sun, 201 1[ 

Daganzo et al., [20 11 

), route 


choices ([Knoop et al.[ |2012[), signal settings and turning movements ([Wu et al.] [201 1[ [Gayah and 


Daganzo[ [201 la[ [Farhi et aT| [201 1[), and loading and unloading processes ([Gayah and Daganzo 


201 lb[ ). Efforts also have been devoted to deriving or approximating the shape of MFDs. Earlier in 
( [Herman and Prigogine[|1979| ), a two-fluid model was introduced to town traffic by splitting vehi¬ 
cles into two groups: moving and stopped vehicles. Relations among flow, density, and speed were 
derived by assuming the average speed in an urban network is a function of the fraction of stopped 


vehicles at any given time. In (Daganzo and Geroliminis| |2008[ ), variational theory was used to 
approximate the MFD for any multi-block, signal-controlled streets without turning movements. 


This method was later extended in (Leclercq and Geroliminis, 2013) to provide an exact solution. 


In (Helbing, 2009), the MFD was derived by averaging the link-based fundamental diagrams for 


all links with the same parameter settings. In ([Daganzo et al.[ [2011[), a two-bin model was used 


to study the possible stationary states in a double-ring network, and multivaluedness was observed 
in the network flow-density relations. However, the aforementioned analytical studies have their 


Daganzo and Geroliminis, 

00 

o 

o 

(N 

Helbing, 2009[), while signal settings were not explicitly modeled 

in ([Herman and Prigogine 

1979 

; Daganzo et al.[ 2011[). Furthermore, the definition of stationary 

states in signalized networ 

■cs has not been clearly defined. 


ized) double-ring network within the framework of kinematic wave models. It was found that 
infinitely many stable states with zero-speed shock waves (ZS) could arise at the same network 
density in the congested branch of MFDs. Simulations using the cell transmission model (CTM) 
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( |Daganzo| |1994| |1995| > of a signalized double ring network unveiled periodic traffic patterns and 
these were defined as stationary states. However, due to the infinite-dimensional state variables 
introduced by the kinematic wave model, traffic statics and dynamics are very difficult to solve 
when traffic signals are present. To the best of our knowledge, the number of stationary states 
in a signalized double-ring network and their stability properties and relations to cycle lengths, 
retaining ratios, and initial densities are all still unclear. 

This study aims to fill this gap by studying a simple signalized intersection shown in Figure 
l(a)[ If exiting vehicles in the downstream are immediately added into the network from their 


corresponding upstream entrances, the signalized intersection is further changed into a signalized 
double-ring network with periodic boundary conditions, which is shown in Figure [T(b)| In the sig¬ 
nalized double-ring network, traffic statics and dynamics are easier to solve because we only need 
to keep track of traffic states on the two rings. Furthermore, instead of using the kinematic wave 



Figure 1: A signalized double-ring network. 


model in ( |Jin et al.[ |2013| >, we use a simplified approximation called the link queue model in ( |Jm 


2012| ) to formulate the traffic dynamics in the signalized double-ring network. With the assump¬ 
tion of a triangular traffic flow fundamental diagram ( |Habermanl|1977| >, the signalized double-ring 
network is studied as a switched affine system. The switching rule is governed by three sets of pa¬ 
rameters: initial densities, route choice behaviors, and signal settings. Periodic density evolution 
orbits are found and defined as stationary states. A Poincare map approach ( [Wiggins et al.H199^ 


Teschl[ |2012| ) is used to study the properties of such stationary states. Poincare maps were origi¬ 
nally used to study the movements of celestial bodies and are a very important tool in analyzing 
periodic orbits: each periodic orbit corresponds to a fixed point on the Poincare map, and its sta¬ 
bility is directly related to the stability of the fixed point. In ( |Jin[ [2013| ), Poincare maps have been 
applied to study the stability and bifurcation in traffic flow in diverge-merge and other networks 
within the framework of network kinematic wave theories. With short cycle lengths, closed-form 
Poincare maps are derived for signalized networks, and stationary states and the corresponding sta¬ 
bilities are obtained by solving and analyzing the fixed points on the Poincare maps. The impacts 
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of retaining ratios and initial densities on the MFDs and the gridlock times are analyzed. With 
long cycle lengths, closed-form Poincare maps are hard to obtain, and thus, the corresponding 
fixed points are numerically solved using the secant method ( [Epperson! |20 1 4[ ) . 

The rest of this paper is organized as follows. In Sectionwe formulate the traffic dynamics 
in a signalized double-ring network as a switched affine system when the link queue model ( |Jin[ 
2012|) and a triangular traffic flow fundamental diagram (jHaberman! |1977|) are used. In Section 


^ we derive the Poincare maps from the density evolution orbits in the switched affine system. 
Stationary states and their stability properties are defined in terms of the fixed points and their 
stability properties on the Poincare maps. In Section we solve the stationary states and the 
corresponding stability properties with short cycle lengths. In Sectionwe discuss the impacts of 
retaining ratios and initial densities on the MFDs and the gridlock times. In Sectionwe provide 
numerical solutions to the fixed points on the Poincare maps when the cycle lengths are long. In 
Section [^ we draw our conclusions with some future research directions. 


2 A switched affine system for a signalized double-ring net¬ 
work 


2.1 A link queue representation 


According to the link queue model (LQM) ( |Jin[ |2012| ), vehicles in a link are treated as always 
being in a queue, and therefore, each link has only one state variable, the average link density. The 
left and the right rings in Figure [T(b^ are denoted as rings 1 and 2 with average densities ki{t) and 
k 2 {t), respectively. We assume both rings have the same length L. Since it is a closed network, 
we have k\(t) +k 2 {t) = 2k, where k is the average network density. In addition, we assume both 
rings have the same location-and-time independent fundamental diagram, qa = Q{ka), a = 1,2, for 
ka e [0,kj], where kj is the jam density. Generally speaking, Q{ka) is a concave function with its 
capacity C obtained at the critical density kc, i.e., C = Q{kc). Then the local demand and supply 
are defined as 


Dait) -- 

= Q{mm{kait),kc}) = 

f Qikait)), 

1 c, 

kait) e [0,kc], 
kait) e ikc,kj\, 

a= 1,2 

(la) 

Sait) -- 

= Qimax{kait),kc}) = 

f c, 

1 Qikait)), 

kait) e [0,kc], 

kait) e ikc,kj], 

a= 1,2 

(lb) 


At the junction, vehicles in the two upstream approaches move alternately, and therefore, there 
are two phases in each cycle. Without loss of generality, we assign phase one to vehicles in ring 1 
and phase two to vehicles in ring 2. The cycle length is T with a lost time A for each phase. The 
green ratio is denoted as K\ for ring 1 and Tt 2 for ring 2. We assume the yellow and all red period 
in each phase is the same as the lost time, and therefore, the effective green time is K\T for ring 1 
and 7t2T for ring 2, and {7ti + 7t2)T = T — 2A. Then the signal regulation can be described by the 
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following two indicator functions: 




1, t G \nT^nT -t-TTjZ'), 

0, otherwise, 

1, t G \nT ^K\T^ {n-\- V'jT — A), 
0, otherwise, 


n G No 

(2a) 

n G No 

(2b) 


where No = {0,1,2,3,According to Equation Q, we have three different combinations of 
(5i(t),^(t)) within one cycle: (i) (1,0) stands for the effective green time period in phase one; 
(ii) (0,1) stands for the effective green time period in phase two; (iii) (0,0) stands for the lost time 
period in either of the phases. 

Due to the signal regulations, the signalized 2x2 junction is equivalent to two alternating 
diverging junctions, and the invariant first-in-first-out (FIFO) diverging model ( |Daganzo[ |1995| > is 
used. The retaining ratio is denoted as (t) G (0,1) for ring 1 and ^2(0 G (0,1) for ring 2, which 
means the turning ratio is 1 — (^i (t) from ring 1 to ring 2 and 1 — (§2(0 from ring 2 to ring 1. Then 
the out-fluxes gi(0, ^2(0’ the in-fluxes /2(0 be calculated as 



(3a) 

M(() = &(») min{D2((), |||, 

(3b) 

/i(0=gi(0^i(0+g2(0(i-^2(0), 

(3c) 

/2(0 = (0(1-^l(O)+52(0^2(0- 

(3d) 


Since it is a closed network, we only need to focus on the traffic dynamics in one of the rings, 
e.g., ring 1. According to the traffic conservation in ring 1, the following equation holds: 


— = j(/i (t) - siiD) =- 1 - SI (t) + — I —«2«). 

With Equations ([T]) to Q, the system equation can be written as 
dki{t) 


(4) 


dt 


= F{ki (0 ; k, 5i (0 , 52 (0 , ^1 (0 > ^2 (t)), 


(5) 


which is closely related to the average network density k, the signal settings (5i (0,52(0), the 
route choice behaviors ((§i(0,<^2(0)- Note that Equation (|5]) is a nonlinear ordinary differential 
equation (ODE) with periodic forces. It is simpler than the kinematic wave model (Equation (3) 
in ( |Jin et al.] |2013| >), which is a partial differential equation (PDE), but still quite challenging to 
solve. 
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2.2 A switched affine system 

For the remainder of this paper, we adopt the following triangular traffic flow fundamental diagram 
(IHabermanl 119771): 


Q{ka) = mm{vfka,Wa{kj-ka)}, a = 1,2, 


( 6 ) 


where v/ and Wa are the free-flow speed and the shock-wave speed, respectively. Then the signal¬ 
ized double-ring network can be further formulated as a switched affine system ( |E1 Aroudi et al.[ 
2007^ [Su^n and Ge||2011[), and the traffic dynamics in Equation Q can be rewritten as 


dk\{t) 

dt 


— Aik\(t )i — 1,2, 


(7) 


where N is the number of possible combinations of {Ai,Bi). The choice of i (or {Ai,Bi)) is governed 
by three sets of parameters: the initial densities (A:i(r),k), the signal settings (5i(t),52(t)), and the 
route choice behaviors (t), (^2(0)- 

Since densities in both rings vary during the effective green times, the {k \, k) space can be 
divided into different regions, in which the coefficients {Ai,Bi) remain the same. According to 
Equations Q, ([^, and (|^, the {ki,k) space can be divided into four regions when (5i(t), 52(t)) = 
(1,0), which are provided in Figure 2(a)[ The meanings of the bold dashed lines in Figure [2^ are 
explained as follows: 


(i) the line k = ^ 


(ii) the line k 


2kc 


^\kj+{\—^\)kc+k\ 


stands for the case when D\ = -r^ < C < 


stands for the case when D 


1 




1-^1 
5i _ 52 


1-^1 

c < 

C < ^1, 


(iii) the line k = stands for the case when <D\ = C; 

(iv) the line k\ = kj — {kj — kc) stands for the case when Di = ||- = C < 

Similarly, the {ki,k) space can be divided into another four regions when (5i(t),^(t)) = (0,1), 
which are provided in Figure |2(b)[ The meanings of the bold dashed lines in Figure |2(b)| are 
provided below: 


(i) the line k ^ ctonHc fr^r* — 


stands for the case when D 2 = ^ = C < 
stands for the case when < £>2 = C; 


(ii) the line k 


{l-2^2)kj+ki 


2(l-&) 


(iii) the line k = y + 2 (i-^ 2 ){k^-k ) stands for the case when D 2 = ^ ^ ^ 

(iv) the line ki =kj —{I — ^ 2 ){kj — kc) stands for the case when D 2 = 


6 





































(a) (5i(r),52(f)) = (l,0) 



(b) (5i(r),52(r)) = (0,l) 

Figure 2: Regions in the {ki,k) space. 
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During the lost time periods (i.e., (5i(t),^(t)) = (0,0)) in phases one and two, densities in both 
rings remain unchanged, and therefore, we consider them as two different regions with the same 
coefficients in Equation Q (i.e., (A,-,5,) = (0,0)). Thus, we have N = 10. The possible values of 
Ai and 5,- and the corresponding conditions are provided in Table [T| 


Table 1: Possible values of A, and and the corresponding conditions 


Region (/) 

Ai 

Bi 

Conditions 

1 

- 7 i 

0 

( 5 i(r),^(/)) - (1,0), Q<h< kc, ’^l<k< 

2 

0 

- 7 \kc 

-{\, 0 ),kc<k\<kj ^i{kj kc), 1 <k<^^ ^‘ 

3 

n 

-Yikj 

( 5i(r),52(/)) - (1,0), ^i{kj kc) <ki <kj, <k< kj-\-^^^ki 

4 

-73 

(N 

1 

1 

(5i(r),d2(f)) = (1,0), 0 < A:i < kj 

[(l-?l)^y-(2-?lW<:i ^lki+{^-^i)kc+ki 2^i-l, , 1,1 ^ u ^ kj+h 

maxi 2 2kc ^ 2 ^ 2^1 4” 2^1 1J ^ ^ — 2 

5 

-74 

2^74 

(gi= (0,1), 0 < fci < kj, 1 < ^ < min{h+fe, 1 > 

6 

0 

74^c 

(5i(r),52(r)) = (0,1), 0 < fci < kj - (1 - ^i^kj-h), '^<k< 

7 

75 

(N 

1 

(5i(r), ^(0) - (0,1), 0 < jfci < kj, max{^ 0-2fe)q+^i} < 

8 

-n 

Yikj 

(di(t), 52 it)) = ( 0 , 1 ), kj - {1 - ^ 2 )ikj -kc) <ki <kj 
ki 1 kc(kj-ki) ^ ^ (l- 2 ^ 2 )kj+ki 

2 + 2(l-&)(q-^,) - - 2(l-fe) 

9 

0 

0 

(5i(r),52(0) = (O5O), and transition from (5i(r),§2(0) = 

10 

0 

0 

{ 5 \{t),& 2 {t)) = (0,0), and transition from (5i(r),^(r)) = (0,1) 


where 71 = 


(l-<^l)^/ ^ _ (l-<^l)v/^c 
L ’ /2 — L^^(^kj-kc) 


73 


L(kj-kc) ’ "4 


L 


, and 75 = 


L^2(kj-kc) ■ 


3 Periodic density evolution orbits and derivation of Poincare 
maps 

3.1 Periodic density evolution orbits 

Due to periodic signal regulations, the switched affine system periodically visits the (^i, k) space in 
Figure]^ and therefore, density evolution orbits are formed. The circulating period is fixed, which 
is the cycle length T. In Figure |3(a)| we provide the density evolution orbit within one cycle in the 
{k\,k) space. And also, since it is a closed network (i.e., k is fixed once it is given), we can map 
the density evolution orbit to the ki axis, which is shown in Figure [3(b^ The evolution process can 
be described as follows: 

(1) At the beginning of the cycle at time t, the initial densities are denoted as (ki{t),k). During 
[t,t + 7t\T),we have (, ^ 2 ) = (1,0). The density k\ decreases as time elapses since vehicles 
in ring 1 can diverge to either of the rings while vehicles in ring 2 have to wait. The density 
evolution may cross multiple regions during this time period, and i can take multiple values 
from 1 to 4 in Table [U 
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(2) At time t-\-%iT, the system switches to region 9. During [t-\-%iT-\-%iT + A), we have 
(5i, ^ 2 ) = (0,0). The densities in both rings remain the same since no vehicles can move to 
the downstream links. 


(3) At time t -\-Tt\T + A, the system switches to phase two. During \t -\-7t\T + A,t + T — A), we 
have (5i,^) = (0,1). The density ki increases as time elapses since vehicles in ring 2 are 
allowed to diverge to either of the rings while vehicles in ring 1 have to wait. Similarly, the 
density evolution may cross multiple regions during this time period, and i can take multiple 
values from 5 to 8 in Table [H 


(4) At time t + r — A, the system switches to region 10. During [t + T — A,t + T), ( 81 , 82 ) = (0,0) 
and the density in ring 1 remains the same. 


(5) At t + r, the system switches back to Step (1) with (, ^ 2 ) = (1,0). The densities become 
(^1 (t + T),k) and serve as new initial densities in the next cycle. 


Based on the above description, it is possible for the signalized double-ring network to have 
periodic density evolution orbits if ki(t + T) = ki(t). Such periodic patterns were observed in 
the simulations and defined as stationary states in ( |Jin et al.[ |2013| ). To calculate the asymptotic 
average network flow-rate, the following equation was used: 


_ gl(0 +g2(0 _ + Is=t-T S2is)ds 


( 8 ) 


When the system reaches a stationary state, q{t) becomes a constant value, which depends on the 
average network density; i.e., q(t)\t^+o^ = q{k). Then the relation between q{k) and k is the MFD 
for the signalized double-ring network. 


3.2 Derivation of Poincare maps 

In the following sections we study the static and dynamic properties associated with the density 
evolution orbits in terms of Poincare maps. According to the discussion in Section [TTI the density 
evolution orbit for each cycle is combined with the following four local maps, Pj, i= 1 , 2 ,3,4: 


ki{t + TtiT)=Piki{t), (9a) 

k\(t 7t\T + A) = P2k\(t Tt\T), (9b) 

k\(t -\-T — A) = Pjiki(t Tt\T + A), (9c) 

ki{t + T)=PAki{t + T-A). (9d) 


If we define the Poincare section as the time when the system first visits the (ki , k) space in Figure 


3(a) in each cycle, the Poincare map can be derived as the composition of the four local mappings. 


ki(t -\-T) = Pk\ (t) = P 4 0 P 3 0 P 2 oPiA:i(t). 


( 10 ) 
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(a) In the {ki , k) space 


0 


k^{t + n{r) 
- 




fci(t)fci(t + r) kj 




kiit + n,T + A) 


-r 




\ 

4^ 


tiCt + r-A) 



(b) In the mapping of ki 

Figure 3: The density evolution orbit within one cycle. 


The Poincare map in Equation ( [TO] ) is well defined and can be analytically derived if initial densi¬ 
ties, retaining ratios, and signal settings are given. Then we can have the following definition of 
stationary states: 

Definition 3.1 (Stationary states) The signalized double-ring network is in a stationary state 
when there exists a fixed point on the Poincare map Pk\{t) satisfying k'l = Pkf 

Due to the existence of noise in transportation networks, it is necessary to understand the 
stability properties of stationary states. Suppose at the beginning of one cycle, there exists a per¬ 
turbation that drives the density in ring 1 from to ki (t) = + e(t), where e(t) is a small error 
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term. Then after n cycles, the density in ring 1 and the error term become 


(lla) 

(llb) 


k\{t-\-nT) = P^k\{t)^ n G No 

e{t + nT)=P’^kx{t)-k\ = {dPfe{t), n G No 


Following ( |La Salle[|1976[[Tbschl[|2012| ), we have the following definition of stability for the fixed 
points on the the Poincare maps: 


Definition 3.2 (Local stability) The fixed point is: (i) Lyapunov stable if for any given j3 > 0, 
there exists an CO < {5 such that \e{t + nT)\ < /3 for any |£(t)| < CO and n G No; (ii) unstable 
otherwise. Furthermore, k\ is asymptotically stable if it is stable and e{t-\-nT) -^0 as n^ +oo. 


4 Stationary states and their stability properties with short cy¬ 
cle lengths 

For simplicity, we consider the following homogeneous settings in the rest of the paper: both rings 
have the same effective green times, i.e., 7r2 = and the same time-independent retaining 

ratios, i.e., = 

4.1 Possible stationary states 

To derive closed-form Poincare maps, we require the cycle length to be short enough so that it 
guarantees the density evolution orbit doesn’t cross multiple regions during each effective green 
time period, i.e., the density evolution orbit stays only within region / G {1,2,3,4} when (5i, 82 ) = 
(1,0), and only in region j G (5,6,7,8} when (5i,^) = (0,1). Under different retaining ratios, 
e.g., > 0.5 and ^ < 0.5, we have 11 different combinations of regions in one cycle, which are 

provided in Table|^ Here, we define the sum of the coefficients (A,, 5,) and {Aj,Bj) as 


X{k\,k) — Ajk\-\-Bj-\-Aik\FBi. (12) 

When the network is not gridlocked, k\ will keep increasing or decreasing if X{k\,k) is always 
greater or smaller than zero, which means it is impossible to have stationary states with the com¬ 
bination of regions {i,j). But if X{k\,k) can take both positive and negative values or is always 
zero, it is possible to have stationary states inside the combination of regions (/, j). In Tablewe 
provide possible values of X {k \, k) and the corresponding conditions for different combinations of 
regions under different retaining ratios. Specifically, when ^ = 0.5, we only have the following 
combinations of regions {ijf. (1,5), (1,6), (1,7), (2,5), (2,6), (4,7), (3,5), and (3,8). The val¬ 
ues of X{k\,k) and the corresponding conditions in regions (1,5), (1,6), (1,7), (2,5), (2,6), and 
(3,5) are the same as those with ^ < 0.5 in TableHowever, the values of X{ki,k) in regions 
(4,7) and (3,8) are always zero with ^ = 0.5. Therefore, we have the following lemma: 
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Table 2: Possible values of X (^i, k) in the 11 combinations of regions 


0.5 < <§ < 1 

0 < <§ < 0.5 

Regions 

X{kuk) 

Condition 

Regions 

X(ki,k) 

Condition 

(1,5) 

>0 

ki <k 

(1,5) 

>0 

ki <k 

= 0 

II 

^2? 

= 0 

II 

^2? 

<0 

ki> k 

<0 

ki > k 

(1,6) 

>0 


(1,6) 

>0 


(1,7) 

>0 

71 < 75 or 

7 i > 75 and ki < ^ 

(1,7) 

>0 


= 0 

7i > 75 and h - ^ 

<0 

7 i > 75 and ki > ^ 

(2,5) 

<0 


(2,5) 

<0 


(2,6) 

= 0 


(2,6) 

= 0 


(2,7) 

<0 


(4,6) 

>0 


(4,7) 

<0 


(4,7) 

>0 


(3,5) 

>0 

72<74and/:i<^(,^Vv 

kr ^ 

(3,5) 

<0 


= 0 

/t- 

72<74andki= ^ 

kr 1 

<0 

72<74and/:i> ^ 

kc 1 

or > 74 

(3,6) 

>0 


(2,8) 

<0 


(3,7) 

>0 

k\> k 

(4,8) 

>0 

ki < k 

= 0 

II 

-2? 

= 0 

II 

-2? 

<0 

k\ <k 

<0 

k\> k 

(3,8) 

>0 


(3,8) 

<0 



Lemma 4.1 (Possible regions having stationary states) When the network is not gridlocked, it is 
only possible for the following combinations of regions {i,j) to have stationary states: 

(i) (1,5), (1,7), (2,6), (3,5), anJ(3,7)/or0.5<(^ <1; 

(ii) (1,5), (2,6), flnr/(4,8)/or0< <0.5; 

(Hi) (1,5), (2,6), (4,7), and (3, S) for ^ =0.5. 

However, when the network is gridlocked with one or both rings jammed, it is only possible for the 
combinations of regions (4,7) and (3,8) to have stationary states with any ^ G (0,1). 

The proof is simple and thus omitted here. Note that when the network is gridlocked with one or 
both rings jammed, the gridlock states are stationary states since they satisfy Definition |3.1[ It is 
easy to identify that these stationary states are inside the combinations of regions (4,7) and (3,8). 
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After identifying the possible combinations of regions having stationary states, we can ana¬ 
lytically derive the closed-form Poincare maps. The corresponding fixed points are the stationary 
states we are interested in. Therefore, we have the following theorem: 

Theorem 4.2 (Possible stationary states) According to Equation Q and Table [7] the Poincare 
maps and the corresponding fixed points in the possible regions having stationary states are de¬ 
rived and provided in Table^ 


Table 3: Poincare maps and fixed points in the possible regions having stationary states 


Regions 

Poincare map Pk\ (t) with 0.5 < (§ < 1 

Fixed points k\ 

(1,5) 

2/t(l-e-ri^3')+yti(f)e-2ri^3' 

2k 

l+e-ri«- 

(1,7) 

- 1) + ki(t)e(i'5-ri)^r 

{kj-2k){e^5’^‘ -1) 

l g(r5-ri)7rr 

(2,6) 

ki (t) 

ki (t) 

(4,7) 

{kj-2k){e(y^-y^^^^ - 1) + ki(t)e(^5-z>)^7' 

1 

(N 

(3,5) 

kj{l- + 2k{l- + h 


(3,7) 

kj{2ey^^^ - - 1) - 2k(el^^^ - 1) + ki (Oe^K^^r 

2k+kj(e'^2’’'‘-l) 

(3,8) 

kj{\- + ki 

kj 

Regions 

Poincare map Pk\{t) with 0 < <§ < 0.5 

Fixed points k^ 

(1,5) 

2k{\-e-y^^^)+h{t)e-^y^^^ 

2k 

l+e-ri,tT 

(2,6) 

h (t) 

kfit) 

(4,7) 

{kj-2k){e^y>-y^^^'^ - 1) 

2k — kj 

(4,8) 

kj{e-^y^^y - 2e-y^^'^ + 1) - 2k{e-^y^^'^ - 
+k\ {t)e~^y^^^ 

kj(\-e-n’^‘)+2ke-n’^‘ 

l+<,-r3«- 

(3,8) 


kj 

Regions 

Poincare map Pki (t) with = 0.5 

Fixed points k\ 

(1,5) 

2k{l -\-ki 

2k 

\+e-n^ 

(2,6) 

k\ {t) 

k\ {t) 

(4,7) 

k\ {t) 

k\ (0 

(3,8) 

k\ {t) 

k\ (0 


4.2 Stability properties of the stationary states 

In this subsection, we analyze the stability properties of the stationary states provided in Table 
particularly the gridlock states. According to Equation (lib), the error term changes to e{t + T) = 
dPe{t) after one cycle. In Tablewe provide the error term e{t-\-T) in different combinations 
of regions under different retaining ratios. According to Definition |3.2[ we have the following 
theorem. 


Theorem 4.3 (Stability properties of the stationary states) When ^ > 0.5, the stationary states 
are: (i) asymptotically stable in the combinations of regions (1,5), (1,7), (4,7), (3,5) and (3,8); 
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Table 4: Changes in the error term after one cycle 


Regions 

e(t + r) with 0.5 <1 

(1,5) 

e{ty^^ 

(1,7) 

e(?)e(r5-ri)^r 

(2,6) 

e(t) 

(4,7) 

e(?)e(r5-%)^r 

(3,5) 

e(t)e(^-^4)7rr 

(3,7) 

e(t)?^ 

(3,8) 

e(t)e^yi-n)^T 

Regions 

e(t + 7 ) with 0 < (^ < 0.5 

(1,5) 

e{t)^^^ 

(2,6) 

e(t) 

(4,7) 


(4,8) 

e(ty=^^ 

(3,8) 


Regions 

e(t-\-T) with ^ = 0.5 

(1,5) 


(2,6) 

e(t) 

(4,7) 

e(t) 

(3,8) 

e(t) 


(ii) Lyapunov stable in the combination of regions (2,6); (Hi) unstable in the combination of re¬ 
gions (3,7). When ^ < 0.5, the stationary states are: (i) asymptotically stable in the combinations 
of regions (1,5) and (4,8); (ii) Lyapunov stable in the combination of regions (2,6); (Hi) unsta¬ 
ble in the combinations of regions (4,7) and (3,8). When ^ = 0.5, the stationary states are: (i) 
asymptotically stable in the combination of regions (1,5); (ii) Lyapunov stable in the combinations 
of regions (2,6), (4,7) and (3,8). 

Proof: From Tablethe error term e(t + T) can be written as e{t + T) = 8 (t)e‘^, where c is the 
coefficient. To obtain the stability properties, we only need to identify the values of the coefficient 
c. The fixed point is unstable if c > 0, while it is stable if c < 0. Furthermore, the fixed point is 
asymptotically stable if c < 0. Since = ^2 = we can have y \ = 74, 72 = 75. Normally, we 
assume kj ~ 5kc, so we have 71 > 75 when ^ > 0.5. In addition, we have 73 > 75 when ^ > 0.5, 
and 73 < 75 when ^ <0.5. Then the derivation of stability is easy and thus omitted here. ■ 
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5 Analysis of macroscopic fundamental diagrams and gridlock 
times 


5.1 Macroscopic fundamental diagrams 

Once the stationary states under different retaining ratios are obtained, the MFD can be analyti¬ 
cally derived. However, the calculation of the average network flow-rate in Equation ([^ involves 
the integral of out-fluxes gi{t) and g 2 {t), which is difficult to solve. Therefore, we use the fol¬ 
lowing approximation when the signalized double-ring network has reached a stationary state, i.e., 
k\{nT) =k\,nE Nq. 

8i82is)ds _ gi{s)ds 

— 2T — 2T 

_ nT{gi{nT)+gi{nT+nT)) _ ^ gi{k^)+gi{k^i{nT+nT)) ^ ^^ 2 ) 

According to Equation (flsj), we have the following MFD when 0.5 < ^ < 1: 


qik;7t,^) 


/ 


< 




JtVfk^ 

ttC, 


TlC 

kC 


kj- 2 k 

kj-k 


^{kj-kc)' 


0 , 


ke [0,kc] 

ke lkc,kj-^{kj-kc)] 

ke {kj-^{kj-kc),kj) 

ke {kj/2,kj] 


(14) 


The proof is provided in The MFD in Equation ( fT?) ) is shown in Eigure |4(a)[ Erom the 
figure, we find that: (i) multivaluedness exists in the network flow-density relation when k e 
((1 (ii) for k > the network can have stationary states with non-zero flow- 

rates, but these stationary states are either Lyapunov stable or unstable; (iii) the gridlock states are 
asymptotically stable. 

Similarly, when 0 < ^ < 0.5, we have the following MED: 


q{k; 7t,^) 


nvfk, ke[0,kc] 

nC, ke{kc,kj-{l-^){kj-kc)] 

’ ke{kj-{l-^){kj-kc),kj) 

0, ke{kj/2,kj] 


(15) 


The proof is similar to and thus omitted here. The MED in Equation ( fTS] ) is provided in Eigure 
4(b)[ Erom the figure, we find that: (i) the signalized double-ring network can maintain higher 


average network flow-rates when the retaining ratios are small; (ii) multivaluedness also exists 
when k > ky/2; (iii) the gridlock states are unstable with small retaining ratios. 

However, when ^ = 0.5, we find that the initial densities in the combinations of regions (2,6), 
(4,7), and (3,8) are all stationary states. Therefore, the MED is very different from those in 
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Equations ([T4|) and (fTS]) and is provided below: 


q{k-,kx,Tl,^) 


TtVfk, 

TtC, 

j^kj—2k-\-ki 


kj-\- 3 kc 


ke [0,kc 

k e {kc, 

ke{‘^^,kj] and 

ki e [max{^^l^^-^^,2k — kj},max{2k—^^‘^^,k}] 


(16) 


The proof is similar to|7]and thus omitted here. The MFD in Equation (16) is provided in Figure 


4(c) 


From the figure, we find that: (i) there are infinitely many stationary states when k > 


kj-\-3kc 


and thus, the corresponding average network flow-rate covers the whole green region shown in the 
figure; (ii) the gridlock states are Lyapunov stable; (iii) there are no unstable stationary states. 


5.2 Gridlock times 


In Theorem 43} we find that stationary states in the combinations of regions (4,7) and (3,8) 
are stable gridlock states when ^ > 0.5. That means starting from any initial conditions in these 
regions, the signalized double-ring network will finally get gridlocked. The gridlock time can be 
analytically calculated from the Poincare maps. For example, in the combination of regions (3,8), 
we have the Poincare map k\{t-\-T) = Pk\(t) = kj (1 —%)^7’)-|-^i(^)g(')^ For an initial 

state (ki(0),k), we have 


kx{T) = 


(17) 


Then after n cycles, we have 


kx{nT) = kj - {kj-kx{(d))e^'^-'^>^^. 


(18) 


Since 72 < 73 , k\ {nT) will converge to kj as n —)■ -|-oo. If we define the gridlock time Tg as the time 
when ki (Tg) ~ (1 — <y)kj, where <7 is very small, we have 


1 


71(73-72) 


ln{ 


kj-ki(0) 

akj 


}■ 


(19) 


In Figure]^ we provide the gridlock patterns with different retaining ratios and initial densities 
in the combination of regions (3,8). In the figure, the horizontal dashed line is the threshold with 
<7 = 0.01 for gridlock conditions. The vertical dashed lines are the calculated gridlock times from 
Equation ( fT^ . The cycle length is 30s with a lost time of 2s for each phase. The updating time 
step At is 0.01s in LQM simulations. The following two trends can be observed from the figure 
and verified by taking the derivatives with respect to ^ and ki(0) in Equation ( [T^ : (i) given the 
same initial densities, the network is harder to get gridlocked with lower retaining ratios, which is 
shown in Figure |5(a)| (ii) given the same retaining ratios, the network will get gridlocked earlier if 
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ring 1 is initially more congested, which is shown in Figure |5(b)[ The gridlock patterns obtained 
from CTM simulations with At = 0.25s are also provided in Figure]^ We find that the two trends 
can also be found in CTM simulations, which indicates that they are the characteristics of the 
signalized double-ring network itself and are not related to the simulation models. Interestingly, 
we also find that with the same initial settings in Figure the network tends to get gridlocked 
earlier in the CTM. 


6 Numerical solutions to the Poincare maps with long cycle 
lengths 

With long cycle lengths, derivation of closed-form Poincare maps is difficult since it is hard to 
track the time instants when the density evolution orbit crosses the boundaries of different regions. 
However, since the system dynamic equation (Equation (|7])) and the coefficients of different regions 
(Table are known, it is possible to provide numerical methods to solve the Poincare maps. 

Here, we denote <I>(A:i) = ki — Pk \. For given average network density k, retaining ratio |, and 
signal settings (5i(t), ^(t)), finding the fixed points k\ is the same as finding the roots in 4>(ki). 
Note that the right-hand side of Equation Q is a piecewise linear function and is continuous at the 
boundaries, and thus, Pk\ as well as <I>(A:i) is continuous. But since the derivative of ‘I>(A:i) is not 
available, the secant method ([Epperson! |2014[) is used and provided below: 




k’^ — k^ 

Av Av 


-1 


‘I>(k^)-‘I>(k^"^) 


( 20 ) 


is the updated density at step n-|-1. k![, <I>(kj) and <I>(kj“^) are the densities and function 

values at steps n and n — \, respectively. 

In Figure]^ we provide the relations between ^(ki) and ki under different average network 
densities when ^ = 0.55, T = 60s, and A = 4s. From the figure, we know that ^>(ki) can have one 
root, multiple roots, or infinite roots under different average network densities. That means starting 
from different kj, we may get different roots in the end using numerical methods. Therefore, a 
brute-force search in k\ is needed to obtain a full map of the stationary states. In addition, when 
4>(ki) has infinite roots (e.g., when k = 65vpm or 78vpm), judgements on the values of 4>(k”) and 
<I>(k”“^) are needed to avoid zero values in the denominator in Equation (20). The algorithm of 
finding stationary states is provided in|^ 

In Figure [T^ we provide the MFD in the signalized double-ring network with ^ = 0.85, T = 
100s, and A = Os using Equation ( |20| ). As a comparison, we also provide the MFDs obtained from 
the variational theory method and CTM simulations ( jJin et al.[ [2013| ) with the same settings in the 
same figure. From the figure, we find that the MFDs obtained from the three methods are different. 
The variational theory method fails to observe the multivaluedness and gridlock phenomena in 
the MFD. Even though such two phenomena are observed using CTM simulations, the unstable 


congested branch with non-zero average network flow-rates is not observed. Using Equation (20) 
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we can get a complete map of the MFD: the multivaluedness and gridlock phenomena as well as 
the unstable congested branch are all observed. The obtained MFD is similar to the one in Figure 


link level in the link queue model, it is not surprised to find that the flow-rates obtained from the 
Poincare map method are systematically lower than or equal to those obtained from the variational 
theory method and CTM simulations. 

In addition, in Figures |7(b)| and |7(c)[ we provide the MFDs obtained from the Poincare map 
method with ^ = 0.3 and ^ = 0.5, respectively. From the two figures, we find that the MFDs with 
long cycle lengths are very similar to the ones with short cycle lengths in Figures [4(b)| and [4{c)| 


4(a) with small cycle lengths. Since detailed traffic dynamics inside a link are aggregated at the 


7 Conclusions 


In this paper, a signalized double-ring network was studied as a switched affine system under the 
framework of the link queue model ( |Jin[ |2012| ) and the assumption of a triangular traffic flow fun¬ 
damental diagram (|Habermanl|1977[). Stationary states with periodic density evolution orbits were 


shown to exist as the switched affine system periodically visits the {k \, k) space in Figure 3(a) 


Poincare maps were introduced to analyze the existence and the property of these stationary states. 
With short cycle lengths, closed-form Poincare maps were obtained. Stationary states and their 
corresponding stability properties were derived and analyzed by finding and analyzing the fixed 
points on the Poincare maps. It was found that under different combinations of signal settings, 
retaining ratios, and initial densities, stationary states can be Lyapunov stable, asymptotically sta¬ 
ble, or unstable. Macroscopic fundamental diagrams were derived based on the stationary states, 
in which the network flow-rate is a function of the network density, route choice behaviors, and 
signal settings. In addition, the derived MFDs are more complete since the multivaluedness and 
gridlock phenomena as well as the unstable congested branch are all observed. Since the system 
will get gridlocked when k>kj/2 and ^ >0.5, the gridlock time was also analyzed under different 
retaining ratios and initial densities. It was found that the network is harder to get gridlocked with 
lower retaining ratios, but it will get gridlocked earlier if one ring is initially more congested. Since 
closed-form Poincare maps are hard to obtain with long cycle lengths, the secant method was used 
to numerically And the flxed points on the Poincare maps. It was found that the obtained stationary 
states and the MFDs are very similar to those with short cycle lengths. 

Different from the kinematic wave approach in ( |Jin et al.| |2013| ) and the two bin model in 
( [Daganzo et al.l|20TT||Gayah and Daganzo[|20lTb| ), the introduction of the link queue model and 
the switched affine system enables us to analytically study the traffic statics and dynamics in the 
double-ring network with general signal settings. The way of using Poincare maps to find the 
stationary states is physically meaningful and mathematically easier to solve. The findings in this 
paper, e.g., stationary states, stability, macroscopic fundamental diagrams, and gridlock patterns, 
are consistent with those in earlier studies ( [Daganzo et al.[ [201 1[ [Gayah and Daganz^ [201 lb[ |Jm 
et al.H201^ . But the solutions in this paper are more complete. 

From our analysis, when k>kj/2 and ^ > 0.5, there exist two types of stationary states in the 
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signalized double-ring network: the asymptotically stable gridlock states and the unstable states 
with more symmetric density distributions and higher average network flow-rates. To improve 
the network performance, we are interested in developing new signal control and route choice 
strategies to adaptively change the effective green times and retaining ratios. In addition, we are 
also interested in applying the analytical framework in this paper to solve the traffic dynamics in 
more general networks, e.g., a signalized grid network, in the future. 


Appendix A. Approximation of average network flow-rates 

The stationary states are provided in Table|^under different retaining ratios. For ^ > 0.5, the pos¬ 
sible combinations of regions having stationary states are: (1,5), (1,7), (2,6), (4,7), (3,5), (3,7), 
and (3, 8 ). For regions (4,7) and (3, 8 ), the stationary states are gridlock states, and therefore, the 
average network average flow-rates are zero. For the rest of the regions, we can approximate the 
average network flow-rate using Equation ( [T^ . 

(1) For regions (1,5), the fixed point is k\ = Starting with ki{nT) = k\, we can get 

k\{nT -\-kT) = ki{nT)e~^'^^^. Since ring 1 is uncongested, gi{ki) = Vfk\. Therefore, the 
average network flow-rate is 

q{k) 7^ =\KVf{k\^ki{nT+ tiT)) 

= = Ttvfk. ( 21 ) 


(2) For regions (1,7), the fixed point is k\ = Starting with k\ (nT) = kj, we can 

get k\{nT -\-kT) = ki{nT)e~'^^^^. Since ring 1 is uncongested, gi(A:i) = Vfki. In addition, 
since T is small, —yiJtT, and (75 — 'Y\)nT are also small. Therefore, the average 

network flow-rate is 


q(k) ^KVf{k\ -\-ki{nT -\- itT)) = ^KVf(kj — 2k) 


TtVf{kj ■ 


> (71-75^ 


TtC 


l-ete-rp'ir 
{kj- 2 k) 
^{kj-kc)-kc ■ 


( 22 ) 


(3) For regions (2, 6 ), the fixed point is kj = ki (t). In this combination of regions, the out-fluxes 
are restricted by the capacity. Therefore, the average network flow-rate is 


,(k) « ^ g>K)+sm’>T + nT)) ^ 


(23) 


(4) Forregions (3,5), the fixed point is kj = ^ ^ \(ri-r 4 )nT —• Starting with (nT) = 
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k\, we cangetfci(nr -\-nT) = kj{\ — -\-k\{nT)e"^^'^. In this combination of regions, 

the out-flux is governed by the supply in ring 1, i.e., gi{ki) = Therefore, 

the average network flow-rate is 


q{k) 


TtC 


— (iy (1 — )-|-i[ (1)) 

2 ^{kj-kc) 


TtC 


2 kj -2 


2k^Y^-kj^Y2 

74-72 


2 ^ {kj-kc 


= TtC 


{kj-2k)^ 
^ {kj-kc) 


= TtC 


{kj— 2 k) 

^{kj-kc)-kc' 


(24) 


(5) For regions (3,7), the fixed point is kj = • Starting with k\{nT) = kj, we can 

getki(nr -|-;rr) = kj{\ — -\-k\{nT)e'^'^^'^. In this combination of regions, the out-flux 

is governed by the supply in ring 1, i.e., gi(ki) = Therefore, the average 

network flow-rate is 


q{k) 


TtC 


2kj-(kj{\-en’^'^)+k\{\+ey2^T^^^^ 

2^{kj-kc) 


„ 2 kj—{kj + 2 k+kj —1)) 

^^{kj-kc) 


= nC 


kj—k 

^ {kj-kc)' 


(25) 


Based on the average network flow-rates calculated above, the macroscopic fundamental di¬ 
agram for the signalized double-ring network with 0.5 < ^ < \ can be easily derived and thus 
omitted here. ■ 


Appendix B. Algorithm of finding stationary states 

Inputs: k, T, A, and <l>(ki) 

Initialization: vector of stationary states, i.e., SS=[]; minimum value of ki, i.e., ki^min = 
max{2k —ky,0}; maximum value of ki, i.e., ki^max = min{2k,k^}; the threshold of ki, i.e., 
eic', searching step, i,e., Ak; maximum number of iterations, i.e., rimax 

For ki = ki^ffiifi : Ak : ki^ffi^K 

Set kj = k\ and calculate 4>(kj) 
lf^{k^^) ==0 

kj is a root, and add it to SS 

Else 

Set k} = Pkj and calculate 4>(A:}) 

For n=l : Umax 
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If \k\ — A:® I < or $(A:}) == 0, add k\ to SS and break 
^tmp ^ ^1 _ ^^^1 ^ [_Az|_]^ and calculate 

= k\ and = <I>(k}) 
k\=k[^Pand^{k\) = ^{k[^P) 

End for 
End if 

End for 
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— ■ Lyapunov stable 



(a) 0.5 < ,* < 1 




(c) ^=0.5 


Figure 4: MFDs for the signalized double-ring network with different retaining ratios. 
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(k^(0),k)=(130,85)vpm 



(a) With different retaining ratios 

^=0.75 



(b) With different initial densities 

Figure 5: The gridlock patterns with different retaining ratios and initial densities. 
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k=20vpm k=65vpm 




k=78vpm k=120vpm 




Figure 6: Relations between and k\ under different average network densities when ^ 

0.55, T = 60s, and A = 4s. 
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^=0.85 &T=100s &A=Os 



k 


(a) ^ = 0.85 


^=0.3&T=100s&A=0s 



^=0.5&T=100s&A=0s 

1000r. 



0 50 100 150 

k 


(C) ^=0.5 


Figure 7: MFDs with T = 100s, A = Os, and different retaining ratios. 
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